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We investigate the zero-temperature and the finite- 
temperature properties of the two-dimensional antiferromag- 
netic quantum spin system composed of the s = 1/2 and s = 1 
spins. The spin excitation spectrum as well as the thermo- 
dynamic quantities are computed by means of the Schwinger- 
boson mean-field theory. We discuss how the magnetic double 
structure with the gapful and the gapless excitations is gener- 
ated, and apply the results systematically to the Haldane gap 
system in a staggered magnetic field and also to the weakly 
coupled ferrimagnetic mixed-spin chains, ft is confirmed that 
the results obtained are consistent with the experiments for 
the quasi-one-dimensional Haldane compounds ifeBaNiOs. 



I. INTRODUCTION 

The Haldane gap systerrfl has been one of the most 
fascinating subjects in condensed matter physics, for 
which extensive experimental and theoretical investiga- 
tions have been providing a variety of new interesting 
phenomena. One of the hot topics is the magnetic 
double structure in the quasi-one-dimensional (ID) Hal- 
dane systeHiSn observed for the rare-earth compounds 
i? 2 BaNi0 5 J3n3 In the case with R = Y 3 +, the system 
has a disordered ground state with the Haldane gap.l If 
R is substituted by other magnetic ions such as Nd 3+ 
and Pr 3+ ,O0 the s = 1/2 spins on the R 3+ sublattice 
order magnetically at low temperatures, thus giving rise 
to the magnetic double structure composed of the Hal- 
dane gap excitation and the gapless excitation induced by 
the long-range order. This system has been theoretically 
treated as the Haldane chains with a static staggered field 
which are induced by—the ordered s = 1/2 spins on the 
R 3+ sublattice so far.E30 In a recent paper ,0 however, it 
has been pointed out that the dynamics of the s = 1/2 
spins may be also important, stimulating us to treat the 
system by a mixed-spin model. Concerning the mag- 
netic double structure, there is another interesting spin 
system of current interest, i.e. the quantum ferrimagi 
netic chain composed of two kinds of mixed spins.cHTtH 
Several compounds have been already, found, which in- 
deed realize the mixed-spin system.E3c3 This has been 
stimulatias-further intensive theoretical studies on this 
subject .Ii-aii3 In particular, it has been pointed out that 
such a ferrimagnetic chain has the double structure for 
the spin excitation spectrum, which controls -the charac- 
teristic properties in the ferrimagnetic chain£j 

The above two subjects, which have been studied ex- 
perimentally in different contexts, should possess the 
common interesting physics, because both systems are 



characterized by the mixture of two kinds of distinct 
spins. Motivated by these hot topics, we here study the 
properties of the antiferromagnetic mixed-spin system in 
detail, for which the s = 1 and s — 1/2 spin chains 
are stacked alternately. We exploit the two-dimensional 
(2D) system as the simplest model which possesses the 
magnetic long-range order at zero temperature. By com- 
puting the dispersion relation and the staggered mag- 
netization by means of the Schwinger-boson mean-field 
theoryjia we clarify how our system generates the mag- 
netic double structure, which naturally interpolates the 
above two interesting spin systems. We also show that 
the results obtained are consistent with the experimental 
findings for quasi-lD Haldane compounds i^BaNiOvBLl 
This paper is organized as follows. We introduce the 
model and then briefly summarize the Schwinger-boson 
techniques in Sec. II. We discuss the zero-temperature 
properties in Sec. Ill, and finally move to the thermody- 
namic properties in Sec. IV. A brief summary is given in 
Sec. V. 



II. SCHWINGER-BOSON MEAN-FIELD THEORY 

Let us consider a 2D mixed-spin model on the square 
lattice, which is described by the following Hamiltonian, 



H = Ji 



2j ' J 2i+rj,2j 



J 2i+n,2j+l 1J 2i,2j+l 



hJ,V 



^ [^S^y • S^2j +JJ + J-i^>2i+l,2j ' ^2i+l,2'j+r) ' (1) 



where 



(Sfj) is the spin operator at the (i, j)-th site 
in the (x, y) plane, and r\ implies the summation to be 
taken over nearest-neighbor sites. All the exchange cou- 
plings Ji, Ji and J3 are assumed to be antiferromagnetic. 
We here focus on the mixed-spin system composed of 
S A = 1/2 and S B — 1 spins, since it is straightforward 
to generalize the results to the arbitrary-spin case. In 
Fig. 0, we have drawn the mixed-spin model schemati- 
cally. The indices A m and B m label the sublattice in the 
unit cell. Note that this mixed-spin model is constructed 
in two ways depending on how we stack the indepen- 
dent spin chains by introducing three kinds of the cou- 
pling constants. By taking J 3 ^> Ji, J 2 , we can study the 
characteristic properties of the s — 1 Haldane spin chains 
coupled with the s = 1/2 gapless spin chains, which may 
have the relevance to the Haldane gap system in a stag=. 
gered field observed in the compounds ^BaNiOs.crtl 
On the other hand, by setting the coupling constants 



1 



A2 / B2. 
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FIG. 1. 2D mixed-spin model on a square lattice. The 
small (large) circle represents S = 1/2 (S B = 1) spin. The 
bold dashed, the thin dashed and the bold solid lines represent 
the coupling constants Ji, J2 and J3, respectively. Ai,A2,Bi 
and B2 specify four spins in the unit cell. 



Ji J2 j J3 j we can investigate how the independent fer- 
rimagnetic chains (J2 = J^rir^i which have been studied 
extensively in recent years,E^rl!-3 are combined to form the 
2D system. The advantage of our mixed-spin approach is 
that we can systematically describe these two interesting 
systems by continuously varying the parameters J\ , J2 
and J3. 

We employ the Schwinger-boson mean-field -theory 
(SBMFT)EJ to study the above mixed -spin model.Ej This 
method was applied successfully to the 2D Heisenberg 
model with uniform spinsJifl'Ej and theft, to the bilayer 
system,u the double-exchange systemJB the ferrimag- 
netic chainjlj etc. It is known that the SBMFT can de- 
scribe the magnetically ordered phase, which is charae-. 
terized by the condensation of the Schwinger-bosonsHa 
In the Schwinger-boson representation, the spin opera- 
tors are expressed in terms of the boson creation and 
annihilation operators 7J, j a , with the Pauli matrix <x 
as S = |7^<r Q( 37/3 (a, (3 =|, J.). Since the unit cell in our 
model includes four sites, we necessarily introduce eight 
kinds of Bose operators 7 Q = a^ 1 ' , a~a\ ba\ which 
belong to the Ai,A 2 ,Bi,B 2 sublattices, respectively. By 
imposing the constraint 7,-i-7i| + 7ji7i4, = %S A or 2S B on 
each site, we can correctly map the original spin system 
to the boson system. Introducing the Lagrange multipli- 



ers A^ and A,^, the Hamiltonian with the constraints is 
recast to 



H 



E 



E ^ EW-w 



+ E ^ E«-^ fl 

+ J2[Ji(s Al s Bl + s A2 s B2 ) 

+ J 2 S Al S A2 + J 3 S Bl S B2 ] , 
where four bond operators are introduced as 



A ■ — n yl> - n (±1 h y 

~ u 2i,2jT u 2i+r),2jl u 2i,2jl u 2i+rj,2jp 



LL 2i,2 

_ (2) 
~~ a 2i,2 

-„« 



B - a y "' 6 (2) -a yz> b y 

J->ijrj — U "2i,2j+l^ u 2i+ri,2j+ll u 2i,2j+l| 2i+?7,2.j + lT ' 
,( 2 ) 



,(1) *,(!) 
J 2i 

-,(2) 



(2) 

(3) 
(4) 
(5) 
(6) 
_i by 

taking the thermal average < Aij n >= A, < Xfj >= Xa, 
etc., which means that these values are assumed to be 
uniform and static. 

By diagonalizing the mean-field Hamiltonian via the 
Bogoliubov transformation, we have 



*n {L) y^j (1) (2) 

Lyij v — a 2i,2jT a 2i,2j+r]l ~ a 2i,2jl a 2i,2j+i]} ' 



Di 



>,(!) 



',(2) 



,C2) 



J 2i+1.2j] u 2i+1.2j+vl 



We perform a Hartree-Fock decomposition of eq. 



tt _ [ „(n)t An) p(2]»(i>)tflW 



n=l 



(7) 



where a and j3 are the Bose operators for normal modes. 
The corresponding energy spectrums read 



Eq — \fE\ 



E, 



(2) 



Eq 



<E~1 



(8) 



where Eq and E\ are given by 



E 
Ei 



ft- 



\ A + \ Z B - 2di -ei- ^ 
{X\ -Xl-ei + flf 
-4dl ((X A - A B ) 2 - (e k + / k ) 



(9) 



with 



/k 



2 A Ji cos k x 
2DJ3 cos k v 



ek = 2CJ2 cos k, 



By minimizing the free energy thus obtained at finite 
temperatures, we end up with the self-consistent equa- 
tions for A — B,C,D,Xa,Xb, 



1 + 2S J 



E 



^cotlr.if^ 



(n) 
k 



7T 



dX A : 

(n) 



2S B = y2[^ cothKE in)dE^ 

^ J it 2 k dX B 



-8J1A 



E 



dk MdE^ 



(10) 

(11) 
(12) 



2 



-4J 2 C 



-4J,D = 



dk 



— coth 

2 K 



O) 



as 



(n) 



7r 2 k <9z? 



(n) 
k 



(13) 
(14) 



We shall solve these equations numerically for given cou- 
pling constants J\, J2 and J3. 



A. Dispersion relation 



with k = 1/(2/cbT), where we have assumed that the 
bond operators which link the s = 1 and s = 1/2 spins 
take the same mean value, A = B. This completes our 
formulation based on the SBMFT. In the following sec- 
tions, we solve these self-consistent equations to estimate 
the excitation spectrum and the thermodynamic quan- 
tities. Since it is not easy to analytically perform the 
Bogoliubov transformation, we numerically diagonalizc 
the mean-field Hamiltonian to compute the energy dis- 
persions E^ and E^ . 



III. PROPERTIES AT ZERO TEMPERATURE 

In order to treat the ground state properties, it should 
be taken into account that the Bogoliubov particles of 
the a branch in eq.(Q) may condense at absolute zero, 
because the excitation energy E^ has its minimal value 

E^' = at k = while the /? bxanch has a finite 
gap even at T = 0. Sarker et a/EJ showed that the 
long-range order is described by the condensation of the 
Schwinger-bosons for the ferromagnetic and antiferro- 
magnetic Heisenberg models. This is also the case for 
our 2D mixed-spin model on a square lattice. 

Suppose that the bosons condense at the states of 
a^|k=o and af^|k=Oi by fictitiously applying an in- 
finitesimal external staggered field to the A and B lat- 
tices. The self-consistent equations at T = , which 
include the Bose condensation, now read, 

r = ^ cothK ^k=oU-oo 



dk 9E< k 



dE^ljdX, 



(15) 



1 + 2S L 
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dd k 
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dev 
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n 

E 



{n) 



dk dE^ 
~^ d\ B 

dk 



(16) 



7T^ 



, 9E ( k n) 
cos k x —- — 



(17) 



dk , dE v 

COS k. 



(») 



dev 



dk , dE { k n) 
— cos k„ 



dfu 



(18) 



(19) 



We start our discussions with the case that may be 
regarded as the Haldane gap system in a staggered field. 
When Jj = 0, the 2D system is completely decoupled into 
the s = 1 massive Haldane chains and the s = 1/2 mass- 
less spin chains. When we introduce the interchain cou- 
plings among them, the ground state has the long-range 
order mainly due to the s = 1/2 spins. The important 
point is that even though we have the long-range order, 
the Haldane- type gapful excitation still exists. There- 
fore, as far as the case with small J\ and J2 is concerned, 
the system is regarded as the one oftgn-.called the Hal- 
dane gap system in a staggered fieldfliZI In Fig. ^, we 
show the dispersion relation calculated for J\ — J2 = 1/2 
and J3 = 1. As seen in this figure, the lower branch 




FIG. 2. The excitation spectrum for the mixed-spin model 
which may be regarded as the Haldane gap in a staggered 
field (Ji = Ji = 1/2, J3 = 1). Note that the Brillouin zone is 
reduced to a quarter of that for the uniform spin case, since 
the unit cell now includes the four sites. 

is gapless with linear dispersion relation, reflecting the 
antifcrromagnetic long-range order, whereas the upper 
optical mode is mainly composed of the Haldane-type 
excitation. It is observed that the dispersion in the k x - 
direction is indeed weak for the optical mode, since the 
interchain coupling is much smaller than the energy gap 
for the optical mode, which confirms that the system in 
Fig. |^ may be regarded as the Haldane gap system in 
a staggered field. Although we have also performed the 
calculation for the cases with smaller J\ and J2, the ob- 
tained dispersion of the optical branch becomes almost 
flat in the fc^-direction, so that we have not shown them 
here. Concerning this limit of small J\, we here make a 
brief comment on the validity of the SBMFT. When J\ 
takes too small value, for which the system is almost de- 
coupled into independent s — 1 and s = 1/2 chains, the 
SBMFT may lead to a pathological result: although the 
Haldane gap for the s — 1 chain is well described by the 
SBMFT, it is not the case for the gapless s = 1/2 chains, 
for which we are left with a gapful phase. Actually, if the 
value of Ji (= J 2) becomes smaller than J1/J3 ~ 10~ 2 , 



3 



we encounter a problem that the present self-consistent 
calculation does not converge, implying that our assump- 
tion for the antiferromagnetic ordered state does not hold 
anymore. Nevertheless, we find that the correct behavior 
with the antiferromagnetic ground state is still obtained 
except for this small parameter region. 

By increasing the coupling parameters J2, J3 continu- 
ously, we naturally enter in the 2D ordered mixed-spin 
system. Note that although during this process the mag- 
netic double structure is kept unchanged in its typical fea- 
ture, the nature of the optical mode is gradually changed 
from the Haldane-gap excitation: i.e. the gapful excita- 
tion may be equally contributed from both the s = 1 and 
s = 1/2 spin sectors. As a reference, we show the disper- 
sion relation for the mixed-spin model with the isotropic 
bonds J\ = J2 = J3 = 1 in Fig. ||. If we further decrease 




FIG. 3. The excitation spectrum for the 2D mixed-spin 
model (Ji = J 2 = J3 = 1). 

the couplings J2 and J3 , the system gradually approaches 
the quasi-lD ferrimagnetic chains with the periodic ar- 
rangement of spins l/2olol/2ol. The dispersion relation 
obtained in the corresponding parameter region is shown 
in Fig. 0. Now the optical mode with a weak disper- 




FIG. 4. The excitation spectrum for the coupled ferrimag- 
netic chains (Ji = 1, J2 = J3 = 1/2). 

sion in the fcy-direction is essentially the same as that 
found for the ferrimagnetic chainJlj In order to clearly 
observe how the gapless dispersion changes its charac- 
ter in the ferrimagnetic-chain limit, we have shown the 
low-energy dispersion relation at k y =0 in Fig. |^. It is 
seen that with the decrease of the coupling constants J2 
and J3, the /c x -linear dependence is gradually changed 
to the k\ dependence characteristic of the ferrimagnetic 
chain except for the small k x region where the 2D an- 
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FIG. 5. The gapless dispersion for the coupled ferrimag- 
netic chains at k y =0. The solid, the dash-dotted and the 
dashed lines correspond to the cases of J2 = J3 = 0.5, 0.1 
and (Ji = 1). 



tiferromagnetic order still gives rise to the fc^-linear de- 
pendence. At J2 = J3 = 0, the system is reduced to the 
isolated ferrimagnetic chains, for which the 2D charac- 
ter completely disappears, and thus the gapless disper- 
sion simply follows the k%, dependence. We note that 
the ferrimagnetic chains have been experimentally real- 
ized in the compounds suqtuas NiCu(C2C>4)2 ■ 4H2O and 



MnCu(pba)(H 2 0) 3 -2H 2 v 
retically-by many groups .Hi 



nd have been studied theo- 
IB Among others, it has been 



reported^ that this system exhibits the dual properties 
consistent with our results: the physical quantities are 
controlled by the s — 1 antiferromagnetic spin chain at 
high temperatures and by the effective s = 1/2 ferromag- 
netic spin chain at low temperatures. 

The above analysis of the excitation spectrum implies 
that-.tjhe quasi-lD Haldane gap system in a staggered 
field0l3_and the quasi-lD weakly coupled ferrimagnetic 
chains Jjjffijl which have been studied in different con- 
texts experimentally, share common interesting physics 
inherent in the mixed-spin systems. In particular, it 
is highly desirable to experimentally observe the mag- 
netic double structure in the excitation spectrum for the 
ferrimagnetic-chain compounds. 



B. Haldane gap in a staggered field 

Let us discuss the case of the Haldane gap system in 
a staggered field in more detail. We here observe how 
the effective staggered field induced by the s = 1/2 
spin chains affects the properties of the Haldane gap, 
by changing the coupling constants J\ and J 2 . For small 
values of J\ and J2, we define the effective staggered field 
on the s = 1 Haldane chain by = J\ < Sf >, where 
< Sf > is the spontaneous staggered magnetization of 
the s = 1/2 spin chain. We numerically estimate the ef- 
fective staggered field, and show the obtained results in 
Fig. H It is seen that the effective staggered field in- 
creases monotonously with the increase of J\ and J2, as 
should be expected. In a similar way, we also compute 
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FIG. 6. The effective staggered field H|t as a function of 
the coupling constant J\ {=J2) where we set J3 = 1. 



the staggered magnetization < S~ > of the s = 1 spin 
chains, and as well as the Haldane gap which is defined 
as the minimum of the excitation energy in the optical 
branch A = E, 



(2) 

k=(f ,o)- 



These quantities are shown in 



Since we are dealing with the 2D model at zero tem- 
perature, our results may not be directly applied to the 
experiments. Nevertheless, we can confirm whether the 
present results are consistent with the-|experiments for 
the rare-garth compounds i^BaNiOsEHj with R = Nd 3+ 
or Pr 3+ .0tj Experimentally, as the temperature is de- 
creased, the system shows the phase transition to the 
magnetically ordered phase. When the temperature is 
further decreased, the staggered moment for the s = 1/2 
sector develops, and thereby gives rise to the increase 
of the staggered field, which, iadeed enhances the mag- 
nitude of the Haldane gapESI'U Also, the staggered mo- 
ment on the Haldane chains increases, as should be ex- 
pected. These characteristic features are consistent with 
our results, and in particular, the qualitative behaviors 
for the staggered magnetization and the Haldane gap 
shown in Figs. JBl and || agree fairly well with experi- 
mental findings .Era 




FIG. 7. Plots of the staggered magnetization Mst as a 
function of the effective staggered magnetic field. 




FIG. 8. Plots of the square of the Haldane gap A as a 
function of the square of the staggered magnetization Mst- 



^ and U rEollowing the way used for the experi- 
analysis,Ba we have plotted the staggered mag- 



Figs. 

mental analys 

netization (the square of the Haldane gap) as a function 
of the effective staggered field (the square of the stag- 
gered magnetization). It is seen that as the staggered 
field is increased, the staggered magnetization as well as 
the Haldane gap A are increased, being consistent with 
the result pointed out for . th e Haldane chain system with 
a static magnetic field. BEJ 



C. Coupled ferrimagnetic chains 

Before closing this section, we briefly discuss the case 
close to the ferrimagnetic chain, which is realized by tak- 
ing the limit of J\ 3> J2 , J3 ■ Even in this one-dimensional 
limit, the system still exhibits the antiferromagnetic or- 
der as far as J2 and J3 take finite values. As mentioned 
above, at J 2 — J3 = 0, the system shows the ferrimag- 
netic order, and hence the dispersion relation .ipt-.thc 
acoustic mode is changed to the quadratic one.EjO In 
this way, the crossover-like behavior is seen in the gapless 
mode, whereas much simpler behavior is observed in the 
gapful mode. In Fig. ^, we display the spin gap A for 
the optical branch in the ferrimagnetic-chain limit. With 
decreasing J2 and J3, the spin gap monotonically de- 
creases, and reaches the value of A=1.778 at J2 = J3 = 0, 
which is very closq-tp that of the quantum Monte Carlo 
method, A=1.767jl3 as already demonstrated by Wu et 
aZ£3 In this way, in the ferrimagnetic-chain limit, the 
SBMFT may provide the reliable estimates of the physi- 
cal quantities even at ouantitative level. As discussed by 
Yamamoto and Fukui,t3 this optical mode is mainly com- 
posed of excitations in the effective s = 1 antiferromag- 
netic spin chain, whereas the acoustic mode is given by 
excitations in the effective s = 1/2 ferromagnetic chain. 



IV. THERMODYNAMIC PROPERTIES AT 
FINITE TEMPERATURES 

We now move to the thermodynamic properties at fi- 
nite temperatures. The self-consistent equations (|io|)- 
( PH ) are solved numerically at finite temperatures to ob- 
tain the thermodynamic quantities. We here briefly sum- 
marize the results obtained. In the previous section, it 
has been shown that even for the set of the parameters 
Jr = J 2 = 1/2, J 3 = 1 (Ji = 1,J 2 = J 3 = 1/2), the 
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FIG. 9. Plots of the spin gap for the optical branch in the 
limit of the coupled ferrimagnetic chains (J2 = J3). 



turn — / J 



(20) 



Xstag — 
i,3'=l,2 

+ 2(-l) i +^ 1 ((S^S^)) q , tl> ] | q ,^ , (21) 

where {{Sf 1 S , ^ 1 )) q , w , etc., are the retarded spin corre- 
lation functions. In Fig. [ll], we plot the effective Curie 



system may be approximately described by the Haldane 
gap system in a staggered field (coupled ferrimagnetic 
chains), so that we shall show the results for these pa- 
rameters below. We start with the effective spin gaps 
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FIG. 10. The spin gaps for the optical and acoustic 
branches as a function of the temperature T: the solid lines 
correspond to the case for J\ — J2 — 1/2, J3 = 1 while the 
dashed lines for Ji = 1, J2 = J3 = 1/2. Note that at higher 
temperatures our SBMF approach breaks down, so that we 
have plotted the data available for each choice of parameters. 

calculated at finite temperatures, which are shown as a 
function of the temperature in Fig. [HJ It is seen that 
both of the two cases exhibit similar temperature de- 
pendence in the spin gap. For each case, there are two 
distinct spin gaps, corresponding respectively to the op- 
tical mode and the acoustic mode, reflecting the double 
structure for the excitation spectrum. Note that both 
of two modes should be massive at Unite temperatures. 
It may be physically sensible to regard these spin gaps 
as the inverse of the correlation lengths. As should be 
expected, the spin gap for the optical mode increases up 
to the zero-temperature value with the decrease of the 
temperature. On the other hand, the spin gap for the 
acoustic mode decreases, leading to the divergent correla- 
tion length which characterizes the magnetically ordered 
state at T = 0. 

The uniform and staggered spin susceptibilities, x un i 
and Xstag, are calculated by using the standard linear- 
response formulae, 
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FIG. 11. Plots of T^stag for the staggered susceptibility 
and Txuni for the uniform susceptibility as a function of the 
temperature T: the solid lines for Ji = J2 = 1/2, J3 = 1 and 
the dashed lines for J\ = 1, J% = J3 = 1/2 

constants Txuni and T\ stag as a function of T. It is 
seen that the effective Curie constant for the uniform 
sector is gradually decreased as the temperature is de- 
creased, while that for the staggered sector diverges, im- 
plying that the antiferromagnetic correlation is enhanced 
at low temperatures. 

As seen from the above results, the thermodynamic 
properties at finite temperatures show quite similar be- 
havior both for the Haldane chain in a staggered field 
and for the quasi-lD ferrimagnetic chain, reflecting the 
magnetic double structure inherent in these mixed-spin 
systems. 



V. SUMMARY 

We have studied the 2D mixed-spin model for which 
the s = 1/2 and s — 1 spin chains are stacked alternately. 
This mixed-spin model includes two interesting spin sys- 
tems addressed recently: the Haldane gap system in a 
staggered field as well as the ferrimagnetic chain. By cal- 
culating the dispersion relation and the thermodynamic 



G 



quantities by means of the Schwinger-boson mean-field 
theory, we have discussed the magnetic double structure 
inherent in our mixed-spin systems. In particular, we 
have treated systematically the quasi-lD Haldane gap 
system in a staggered magnetic field and also the mixed- 
spin chain with the ferrimagnetic ground state. This 
implies that these two spin systems, which have been 
studied in different contexts experimentally, should pos- 
sess interesting physics common to the mixed spin sys- 
tems. We have also found that the results obtained for 
the staggered- field effect on the Haldane gap system are 
qualitatively consistent with the experimental findings 
in the rare-earth compounds -fi^BaNiOs. It remains an 
interesting problem to evaluate dynamical quantities re- 
lated to the neutron scattering, etc. Also, it is impc*-. 
tant to study how the string-order parameter behaves,E2l 
when the system changes from the Haldane system to 
the ferrimagnetic chain. These problems are now under 
consideration. 
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We investigate the zero-temperature and the finite- 
temperature properties of the two-dimensional antiferromag- 
netic quantum spin system composed of the s = 1/2 and s = 1 
spins. The spin excitation spectrum as well as the thermo- 
dynamic quantities are computed by means of the Schwinger- 
boson mean-field theory. We discuss how the magnetic double 
structure with the gapful and the gapless excitations is gener- 
ated, and apply the results systematically to the Haldane gap 
system in a staggered magnetic field and also to the weakly 
coupled ferrimagnetic mixed-spin chains, ft is confirmed that 
the results obtained are consistent with the experiments for 
the quasi-one-dimensional Haldane compounds ifeBaNiOs. 



I. INTRODUCTION 

The Haldane gap systerrfl has been one of the most 
fascinating subjects in condensed matter physics, for 
which extensive experimental and theoretical investiga- 
tions have been providing a variety of new interesting 
phenomena. One of the hot topics is the magnetic 
double structure in the quasi-one-dimensional (ID) Hal- 
dane systeHiSn observed for the rare-earth compounds 
i? 2 BaNi0 5 J3n3 In the case with R = Y 3 +, the system 
has a disordered ground state with the Haldane gap.l If 
R is substituted by other magnetic ions such as Nd 3+ 
and Pr 3+ ,O0 the s = 1/2 spins on the R 3+ sublattice 
order magnetically at low temperatures, thus giving rise 
to the magnetic double structure composed of the Hal- 
dane gap excitation and the gapless excitation induced by 
the long-range order. This system has been theoretically 
treated as the Haldane chains with a static staggered field 
which are induced by—the ordered s = 1/2 spins on the 
R 3+ sublattice so far.E30 In a recent paper ,0 however, it 
has been pointed out that the dynamics of the s = 1/2 
spins may be also important, stimulating us to treat the 
system by a mixed-spin model. Concerning the mag- 
netic double structure, there is another interesting spin 
system of current interest, i.e. the quantum ferrimagi 
netic chain composed of two kinds of mixed spins.cHTtH 
Several compounds have been already, found, which in- 
deed realize the mixed-spin system.E3c3 This has been 
stimulatias-further intensive theoretical studies on this 
subject .Ii-aii3 In particular, it has been pointed out that 
such a ferrimagnetic chain has the double structure for 
the spin excitation spectrum, which controls -the charac- 
teristic properties in the ferrimagnetic chain£j 

The above two subjects, which have been studied ex- 
perimentally in different contexts, should possess the 
common interesting physics, because both systems are 



characterized by the mixture of two kinds of distinct 
spins. Motivated by these hot topics, we here study the 
properties of the antiferromagnetic mixed-spin system in 
detail, for which the s = 1 and s — 1/2 spin chains 
are stacked alternately. We exploit the two-dimensional 
(2D) system as the simplest model which possesses the 
magnetic long-range order at zero temperature. By com- 
puting the dispersion relation and the staggered mag- 
netization by means of the Schwinger-boson mean-field 
theoryjia we clarify how our system generates the mag- 
netic double structure, which naturally interpolates the 
above two interesting spin systems. We also show that 
the results obtained are consistent with the experimental 
findings for quasi-lD Haldane compounds i^BaNiOvBLl 
This paper is organized as follows. We introduce the 
model and then briefly summarize the Schwinger-boson 
techniques in Sec. II. We discuss the zero-temperature 
properties in Sec. Ill, and finally move to the thermody- 
namic properties in Sec. IV. A brief summary is given in 
Sec. V. 



II. SCHWINGER-BOSON MEAN-FIELD THEORY 

Let us consider a 2D mixed-spin model on the square 
lattice, which is described by the following Hamiltonian, 



H = Ji 



2j ' J 2i+rj,2j 



J 2i+n,2j+l 1J 2i,2j+l 



hJ,V 



^ [^S^y • S^2j +JJ + J-i^>2i+l,2j ' ^2i+l,2'j+r) ' (1) 



where 



(Sfj) is the spin operator at the (i, j)-th site 
in the (x, y) plane, and r\ implies the summation to be 
taken over nearest-neighbor sites. All the exchange cou- 
plings Ji, Ji and J3 are assumed to be antiferromagnetic. 
We here focus on the mixed-spin system composed of 
S A = 1/2 and S B — 1 spins, since it is straightforward 
to generalize the results to the arbitrary-spin case. In 
Fig. 0, we have drawn the mixed-spin model schemati- 
cally. The indices A m and B m label the sublattice in the 
unit cell. Note that this mixed-spin model is constructed 
in two ways depending on how we stack the indepen- 
dent spin chains by introducing three kinds of the cou- 
pling constants. By taking J 3 ^> Ji, J 2 , we can study the 
characteristic properties of the s — 1 Haldane spin chains 
coupled with the s = 1/2 gapless spin chains, which may 
have the relevance to the Haldane gap system in a stag=. 
gered field observed in the compounds ^BaNiOs.crtl 
On the other hand, by setting the coupling constants 



1 



A2 / B2. 




- - Jl J2 J3 

FIG. 1. 2D mixed-spin model on a square lattice. The 
small (large) circle represents S = 1/2 (S B = 1) spin. The 
bold dashed, the thin dashed and the bold solid lines represent 
the coupling constants Ji, J2 and J3, respectively. Ai,A2,Bi 
and B2 specify four spins in the unit cell. 



Ji J2 j J3 j we can investigate how the independent fer- 
rimagnetic chains (J2 = J^rir^i which have been studied 
extensively in recent years,E^rl!-3 are combined to form the 
2D system. The advantage of our mixed-spin approach is 
that we can systematically describe these two interesting 
systems by continuously varying the parameters J\ , J2 
and J3. 

We employ the Schwinger-boson mean-field -theory 
(SBMFT)EJ to study the above mixed -spin model.Ej This 
method was applied successfully to the 2D Heisenberg 
model with uniform spinsJifl'Ej and theft, to the bilayer 
system,u the double-exchange systemJB the ferrimag- 
netic chainjlj etc. It is known that the SBMFT can de- 
scribe the magnetically ordered phase, which is charae-. 
terized by the condensation of the Schwinger-bosonsHa 
In the Schwinger-boson representation, the spin opera- 
tors are expressed in terms of the boson creation and 
annihilation operators 7J, j a , with the Pauli matrix <x 
as S = |7^<r Q( 37/3 (a, (3 =|, J.). Since the unit cell in our 
model includes four sites, we necessarily introduce eight 
kinds of Bose operators 7 Q = a^ 1 ' , a~a\ ba\ which 
belong to the Ai,A 2 ,Bi,B 2 sublattices, respectively. By 
imposing the constraint 7,-i-7i| + 7ji7i4, = %S A or 2S B on 
each site, we can correctly map the original spin system 
to the boson system. Introducing the Lagrange multipli- 



ers A^ and A,^, the Hamiltonian with the constraints is 
recast to 



H 



E 



E ^ EW-w 



+ E ^ E«-^ fl 

+ J2[Ji(s Al s Bl + s A2 s B2 ) 

+ J 2 S Al S A2 + J 3 S Bl S B2 ] , 
where four bond operators are introduced as 



A ■ — n yl> - n (±1 h y 

~ u 2i,2jT u 2i+r),2jl u 2i,2jl u 2i+rj,2jp 



LL 2i,2 

_ (2) 
~~ a 2i,2 

-„« 



B - a y "' 6 (2) -a yz> b y 

J->ijrj — U "2i,2j+l^ u 2i+ri,2j+ll u 2i,2j+l| 2i+?7,2.j + lT ' 
,( 2 ) 



,(1) *,(!) 
J 2i 

-,(2) 



(2) 

(3) 
(4) 
(5) 
(6) 
_i by 

taking the thermal average < Aij n >= A, < Xfj >= Xa, 
etc., which means that these values are assumed to be 
uniform and static. 

By diagonalizing the mean-field Hamiltonian via the 
Bogoliubov transformation, we have 



*n {L) y^j (1) (2) 

Lyij v — a 2i,2jT a 2i,2j+r]l ~ a 2i,2jl a 2i,2j+i]} ' 



Di 



>,(!) 



',(2) 



,C2) 



J 2i+1.2j] u 2i+1.2j+vl 



We perform a Hartree-Fock decomposition of eq. 



tt _ [ „(n)t An) p(2]»(i>)tflW 



n=l 



(7) 



where a and j3 are the Bose operators for normal modes. 
The corresponding energy spectrums read 



Eq — \fE\ 



E, 



(2) 



Eq 



<E~1 



(8) 



where Eq and E\ are given by 



E 
Ei 



ft- 



\ A + \ Z B - 2di -ei- ^ 
{X\ -Xl-ei + flf 
-4dl ((X A - A B ) 2 - (e k + / k ) 



(9) 



with 



/k 



2 A Ji cos k x 
2DJ3 cos k v 



ek = 2CJ2 cos k, 



By minimizing the free energy thus obtained at finite 
temperatures, we end up with the self-consistent equa- 
tions for A — B,C,D,Xa,Xb, 



1 + 2S J 



E 



^cotlr.if^ 



(n) 
k 



7T 



dX A : 

(n) 



2S B = y2[^ cothKE in)dE^ 

^ J it 2 k dX B 



-8J1A 



E 



dk MdE^ 



(10) 

(11) 
(12) 



2 



-4J 2 C 



-4J,D = 



dk 



— coth 

2 K 



O) 



as 



(n) 



7r 2 k <9z? 



(n) 
k 



(13) 
(14) 



We shall solve these equations numerically for given cou- 
pling constants J\, J2 and J3. 



A. Dispersion relation 



with k = 1/(2/cbT), where we have assumed that the 
bond operators which link the s = 1 and s = 1/2 spins 
take the same mean value, A = B. This completes our 
formulation based on the SBMFT. In the following sec- 
tions, we solve these self-consistent equations to estimate 
the excitation spectrum and the thermodynamic quan- 
tities. Since it is not easy to analytically perform the 
Bogoliubov transformation, we numerically diagonalizc 
the mean-field Hamiltonian to compute the energy dis- 
persions E^ and E^ . 



III. PROPERTIES AT ZERO TEMPERATURE 

In order to treat the ground state properties, it should 
be taken into account that the Bogoliubov particles of 
the a branch in eq.(Q) may condense at absolute zero, 
because the excitation energy E^ has its minimal value 

E^' = at k = while the /? bxanch has a finite 
gap even at T = 0. Sarker et a/EJ showed that the 
long-range order is described by the condensation of the 
Schwinger-bosons for the ferromagnetic and antiferro- 
magnetic Heisenberg models. This is also the case for 
our 2D mixed-spin model on a square lattice. 

Suppose that the bosons condense at the states of 
a^|k=o and af^|k=Oi by fictitiously applying an in- 
finitesimal external staggered field to the A and B lat- 
tices. The self-consistent equations at T = , which 
include the Bose condensation, now read, 

r = ^ cothK ^k=oU-oo 



dk 9E< k 



dE^ljdX, 



(15) 



1 + 2S L 



BE*? 



-4A = r 



-2C = r 



-2D = r 



8EV 



dd k 



8R 



dev 



8EU 



k=0 



k=0 



k=0 



k=0 



E 

n 

E 

n 

E 

n 

E 



{n) 



dk dE^ 
~^ d\ B 

dk 



(16) 



7T^ 



, 9E ( k n) 
cos k x —- — 



(17) 



dk , dE v 

COS k. 



(») 



dev 



dk , dE { k n) 
— cos k„ 



dfu 



(18) 



(19) 



We start our discussions with the case that may be 
regarded as the Haldane gap system in a staggered field. 
When Jj = 0, the 2D system is completely decoupled into 
the s = 1 massive Haldane chains and the s = 1/2 mass- 
less spin chains. When we introduce the interchain cou- 
plings among them, the ground state has the long-range 
order mainly due to the s = 1/2 spins. The important 
point is that even though we have the long-range order, 
the Haldane- type gapful excitation still exists. There- 
fore, as far as the case with small J\ and J2 is concerned, 
the system is regarded as the one oftgn-.called the Hal- 
dane gap system in a staggered fieldfliZI In Fig. ^, we 
show the dispersion relation calculated for J\ — J2 = 1/2 
and J3 = 1. As seen in this figure, the lower branch 




FIG. 2. The excitation spectrum for the mixed-spin model 
which may be regarded as the Haldane gap in a staggered 
field (Ji = Ji = 1/2, J3 = 1). Note that the Brillouin zone is 
reduced to a quarter of that for the uniform spin case, since 
the unit cell now includes the four sites. 

is gapless with linear dispersion relation, reflecting the 
antifcrromagnetic long-range order, whereas the upper 
optical mode is mainly composed of the Haldane-type 
excitation. It is observed that the dispersion in the k x - 
direction is indeed weak for the optical mode, since the 
interchain coupling is much smaller than the energy gap 
for the optical mode, which confirms that the system in 
Fig. |^ may be regarded as the Haldane gap system in 
a staggered field. Although we have also performed the 
calculation for the cases with smaller J\ and J2, the ob- 
tained dispersion of the optical branch becomes almost 
flat in the fc^-direction, so that we have not shown them 
here. Concerning this limit of small J\, we here make a 
brief comment on the validity of the SBMFT. When J\ 
takes too small value, for which the system is almost de- 
coupled into independent s — 1 and s = 1/2 chains, the 
SBMFT may lead to a pathological result: although the 
Haldane gap for the s — 1 chain is well described by the 
SBMFT, it is not the case for the gapless s = 1/2 chains, 
for which we are left with a gapful phase. Actually, if the 
value of Ji (= J 2) becomes smaller than J1/J3 ~ 10~ 2 , 



3 



we encounter a problem that the present self-consistent 
calculation does not converge, implying that our assump- 
tion for the antiferromagnetic ordered state does not hold 
anymore. Nevertheless, we find that the correct behavior 
with the antiferromagnetic ground state is still obtained 
except for this small parameter region. 

By increasing the coupling parameters J2, J3 continu- 
ously, we naturally enter in the 2D ordered mixed-spin 
system. Note that although during this process the mag- 
netic double structure is kept unchanged in its typical fea- 
ture, the nature of the optical mode is gradually changed 
from the Haldane-gap excitation: i.e. the gapful excita- 
tion may be equally contributed from both the s = 1 and 
s = 1/2 spin sectors. As a reference, we show the disper- 
sion relation for the mixed-spin model with the isotropic 
bonds J\ = J2 = J3 = 1 in Fig. ||. If we further decrease 




FIG. 3. The excitation spectrum for the 2D mixed-spin 
model (Ji = J 2 = J3 = 1). 

the couplings J2 and J3 , the system gradually approaches 
the quasi-lD ferrimagnetic chains with the periodic ar- 
rangement of spins l/2olol/2ol. The dispersion relation 
obtained in the corresponding parameter region is shown 
in Fig. 0. Now the optical mode with a weak disper- 




FIG. 4. The excitation spectrum for the coupled ferrimag- 
netic chains (Ji = 1, J2 = J3 = 1/2). 

sion in the fcy-direction is essentially the same as that 
found for the ferrimagnetic chainJlj In order to clearly 
observe how the gapless dispersion changes its charac- 
ter in the ferrimagnetic-chain limit, we have shown the 
low-energy dispersion relation at k y =0 in Fig. |^. It is 
seen that with the decrease of the coupling constants J2 
and J3, the /c x -linear dependence is gradually changed 
to the k\ dependence characteristic of the ferrimagnetic 
chain except for the small k x region where the 2D an- 
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FIG. 5. The gapless dispersion for the coupled ferrimag- 
netic chains at k y =0. The solid, the dash-dotted and the 
dashed lines correspond to the cases of J2 = J3 = 0.5, 0.1 
and (Ji = 1). 



tiferromagnetic order still gives rise to the fc^-linear de- 
pendence. At J2 = J3 = 0, the system is reduced to the 
isolated ferrimagnetic chains, for which the 2D charac- 
ter completely disappears, and thus the gapless disper- 
sion simply follows the k%, dependence. We note that 
the ferrimagnetic chains have been experimentally real- 
ized in the compounds suqtuas NiCu(C2C>4)2 ■ 4H2O and 



MnCu(pba)(H 2 0) 3 -2H 2 v 
retically-by many groups .Hi 



nd have been studied theo- 
IB Among others, it has been 



reported^ that this system exhibits the dual properties 
consistent with our results: the physical quantities are 
controlled by the s — 1 antiferromagnetic spin chain at 
high temperatures and by the effective s = 1/2 ferromag- 
netic spin chain at low temperatures. 

The above analysis of the excitation spectrum implies 
that-.tjhe quasi-lD Haldane gap system in a staggered 
field0l3_and the quasi-lD weakly coupled ferrimagnetic 
chains Jjjffijl which have been studied in different con- 
texts experimentally, share common interesting physics 
inherent in the mixed-spin systems. In particular, it 
is highly desirable to experimentally observe the mag- 
netic double structure in the excitation spectrum for the 
ferrimagnetic-chain compounds. 



B. Haldane gap in a staggered field 

Let us discuss the case of the Haldane gap system in 
a staggered field in more detail. We here observe how 
the effective staggered field induced by the s = 1/2 
spin chains affects the properties of the Haldane gap, 
by changing the coupling constants J\ and J 2 . For small 
values of J\ and J2, we define the effective staggered field 
on the s = 1 Haldane chain by = J\ < Sf >, where 
< Sf > is the spontaneous staggered magnetization of 
the s = 1/2 spin chain. We numerically estimate the ef- 
fective staggered field, and show the obtained results in 
Fig. H It is seen that the effective staggered field in- 
creases monotonously with the increase of J\ and J2, as 
should be expected. In a similar way, we also compute 
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FIG. 6. The effective staggered field H|t as a function of 
the coupling constant J\ {=J2) where we set J3 = 1. 



the staggered magnetization < S~ > of the s = 1 spin 
chains, and as well as the Haldane gap which is defined 
as the minimum of the excitation energy in the optical 
branch A = E, 



(2) 

k=(f ,o)- 



These quantities are shown in 



Since we are dealing with the 2D model at zero tem- 
perature, our results may not be directly applied to the 
experiments. Nevertheless, we can confirm whether the 
present results are consistent with the-|experiments for 
the rare-garth compounds i^BaNiOsEHj with R = Nd 3+ 
or Pr 3+ .0tj Experimentally, as the temperature is de- 
creased, the system shows the phase transition to the 
magnetically ordered phase. When the temperature is 
further decreased, the staggered moment for the s = 1/2 
sector develops, and thereby gives rise to the increase 
of the staggered field, which, iadeed enhances the mag- 
nitude of the Haldane gapESI'U Also, the staggered mo- 
ment on the Haldane chains increases, as should be ex- 
pected. These characteristic features are consistent with 
our results, and in particular, the qualitative behaviors 
for the staggered magnetization and the Haldane gap 
shown in Figs. JBl and || agree fairly well with experi- 
mental findings .Era 




FIG. 7. Plots of the staggered magnetization Mst as a 
function of the effective staggered magnetic field. 




FIG. 8. Plots of the square of the Haldane gap A as a 
function of the square of the staggered magnetization Mst- 



^ and U rEollowing the way used for the experi- 
analysis,Ba we have plotted the staggered mag- 



Figs. 

mental analys 

netization (the square of the Haldane gap) as a function 
of the effective staggered field (the square of the stag- 
gered magnetization). It is seen that as the staggered 
field is increased, the staggered magnetization as well as 
the Haldane gap A are increased, being consistent with 
the result pointed out for . th e Haldane chain system with 
a static magnetic field. BEJ 



C. Coupled ferrimagnetic chains 

Before closing this section, we briefly discuss the case 
close to the ferrimagnetic chain, which is realized by tak- 
ing the limit of J\ 3> J2 , J3 ■ Even in this one-dimensional 
limit, the system still exhibits the antiferromagnetic or- 
der as far as J2 and J3 take finite values. As mentioned 
above, at J 2 — J3 = 0, the system shows the ferrimag- 
netic order, and hence the dispersion relation .ipt-.thc 
acoustic mode is changed to the quadratic one.EjO In 
this way, the crossover-like behavior is seen in the gapless 
mode, whereas much simpler behavior is observed in the 
gapful mode. In Fig. ^, we display the spin gap A for 
the optical branch in the ferrimagnetic-chain limit. With 
decreasing J2 and J3, the spin gap monotonically de- 
creases, and reaches the value of A=1.778 at J2 = J3 = 0, 
which is very closq-tp that of the quantum Monte Carlo 
method, A=1.767jl3 as already demonstrated by Wu et 
aZ£3 In this way, in the ferrimagnetic-chain limit, the 
SBMFT may provide the reliable estimates of the physi- 
cal quantities even at ouantitative level. As discussed by 
Yamamoto and Fukui,t3 this optical mode is mainly com- 
posed of excitations in the effective s = 1 antiferromag- 
netic spin chain, whereas the acoustic mode is given by 
excitations in the effective s = 1/2 ferromagnetic chain. 



IV. THERMODYNAMIC PROPERTIES AT 
FINITE TEMPERATURES 

We now move to the thermodynamic properties at fi- 
nite temperatures. The self-consistent equations (|io|)- 
( PH ) are solved numerically at finite temperatures to ob- 
tain the thermodynamic quantities. We here briefly sum- 
marize the results obtained. In the previous section, it 
has been shown that even for the set of the parameters 
Jr = J 2 = 1/2, J 3 = 1 (Ji = 1,J 2 = J 3 = 1/2), the 
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FIG. 9. Plots of the spin gap for the optical branch in the 
limit of the coupled ferrimagnetic chains (J2 = J3). 
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(20) 



Xstag — 
i,3'=l,2 

+ 2(-l) i +^ 1 ((S^S^)) q , tl> ] | q ,^ , (21) 

where {{Sf 1 S , ^ 1 )) q , w , etc., are the retarded spin corre- 
lation functions. In Fig. [ll], we plot the effective Curie 



system may be approximately described by the Haldane 
gap system in a staggered field (coupled ferrimagnetic 
chains), so that we shall show the results for these pa- 
rameters below. We start with the effective spin gaps 
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FIG. 10. The spin gaps for the optical and acoustic 
branches as a function of the temperature T: the solid lines 
correspond to the case for J\ — J2 — 1/2, J3 = 1 while the 
dashed lines for Ji = 1, J2 = J3 = 1/2. Note that at higher 
temperatures our SBMF approach breaks down, so that we 
have plotted the data available for each choice of parameters. 

calculated at finite temperatures, which are shown as a 
function of the temperature in Fig. [HJ It is seen that 
both of the two cases exhibit similar temperature de- 
pendence in the spin gap. For each case, there are two 
distinct spin gaps, corresponding respectively to the op- 
tical mode and the acoustic mode, reflecting the double 
structure for the excitation spectrum. Note that both 
of two modes should be massive at Unite temperatures. 
It may be physically sensible to regard these spin gaps 
as the inverse of the correlation lengths. As should be 
expected, the spin gap for the optical mode increases up 
to the zero-temperature value with the decrease of the 
temperature. On the other hand, the spin gap for the 
acoustic mode decreases, leading to the divergent correla- 
tion length which characterizes the magnetically ordered 
state at T = 0. 

The uniform and staggered spin susceptibilities, x un i 
and Xstag, are calculated by using the standard linear- 
response formulae, 
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FIG. 11. Plots of T^stag for the staggered susceptibility 
and Txuni for the uniform susceptibility as a function of the 
temperature T: the solid lines for Ji = J2 = 1/2, J3 = 1 and 
the dashed lines for J\ = 1, J% = J3 = 1/2 

constants Txuni and T\ stag as a function of T. It is 
seen that the effective Curie constant for the uniform 
sector is gradually decreased as the temperature is de- 
creased, while that for the staggered sector diverges, im- 
plying that the antiferromagnetic correlation is enhanced 
at low temperatures. 

As seen from the above results, the thermodynamic 
properties at finite temperatures show quite similar be- 
havior both for the Haldane chain in a staggered field 
and for the quasi-lD ferrimagnetic chain, reflecting the 
magnetic double structure inherent in these mixed-spin 
systems. 



V. SUMMARY 

We have studied the 2D mixed-spin model for which 
the s = 1/2 and s — 1 spin chains are stacked alternately. 
This mixed-spin model includes two interesting spin sys- 
tems addressed recently: the Haldane gap system in a 
staggered field as well as the ferrimagnetic chain. By cal- 
culating the dispersion relation and the thermodynamic 



G 



quantities by means of the Schwinger-boson mean-field 
theory, we have discussed the magnetic double structure 
inherent in our mixed-spin systems. In particular, we 
have treated systematically the quasi-lD Haldane gap 
system in a staggered magnetic field and also the mixed- 
spin chain with the ferrimagnetic ground state. This 
implies that these two spin systems, which have been 
studied in different contexts experimentally, should pos- 
sess interesting physics common to the mixed spin sys- 
tems. We have also found that the results obtained for 
the staggered- field effect on the Haldane gap system are 
qualitatively consistent with the experimental findings 
in the rare-earth compounds -fi^BaNiOs. It remains an 
interesting problem to evaluate dynamical quantities re- 
lated to the neutron scattering, etc. Also, it is impc*-. 
tant to study how the string-order parameter behaves,E2l 
when the system changes from the Haldane system to 
the ferrimagnetic chain. These problems are now under 
consideration. 
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